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SUMMARY 



Experiments were conducted with two Delta wings and a streamlined 
body with lifting surfaces to investigate the evolution of turbulent 
trailing vortices in stratified and unstratified water. The vortex 
trajectories were determined as a function of the relative depth D/b^, 
normalized time V^t/b^, and the stratification parameter Nb^/V^, The 
results have shown that the vortices rise only to a finite height as 
they gradually demise under the influence of turbulence, sinusoidal 
instability, and the vortex breakdown. The effect of stratification 
is to reduce the lifespan of vortices and the maximum height attained 
by them. Various approximate analyses have been examined and new ones 
have been proposed to determine the lifespan of the vortices. Finally, 
an exploratory numerical analysis of the two-dimensional, unsteady, 
laminar vortices has been carried out with encouraging results. 
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I. INTRODUCTION 



Trailing vortices formed by a lifting surface traveling through a 
stratified or unstratified fluid medium have been the subject of recent 
interest primarily because of their importance in the determination of 
safe aircraft separation distances for terminal areas. In recent years 
this interest has extended to the trailing vortices generated by the 
control planes of submerged bodies such as a submarine partly because of 
the effect of these vortices on the motion of the submerged body itself 
and partly because of the interaction of these vortices with the surround- 
ing fluid medium. 

The trailing vortices may be generated by any number of lifting sur- 
faces such as an airplane wing, a missile, single or double Delta wings, 
or the control surfaces of a submerged body (sail planes, rudder, and 
stern). The initial organization or roll-up of a trailing vortex from 
the vorticity generated by the lifting surface is very rapid and is 
generally complete within the space of several lifting-surface spans. 

Whereas a single vortex is extremely stable, one or more vortices in 
close proximity exhibit a number of complex interactions and give rise to 
new phenomena which tend'to accelerate the demise of the vortices in- 
volved. When a lifting surface starts from rest, a starting-vortex or 
a Horseshoe vortex is formed since a vortex must be connected either to 
itself or to a rigid boundary. A Horseshoe vortex is also formed when a 
lifting surface accelerates or decelerates rapidly, or undergoes a radical 
change in the angle of attack. When a control surface (e.g., a Delta 
wing) starts from rest and moves with a constant velocity, the Horseshoe 
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vortex closes on itself, breaks away from the rest of the trailing vor- 
tex, links rapidly, and forms a vortex ring. While the first vortex 
ring undergoes highly complex deformations, a second vortex ring is 
formed from the Horseshoe vortex connecting the tail end of the trailing 
vortex. This process does not continue forever for a number of reasons. 

The speed of the generating surface is much greater than the speed of 
advance of the vortex rings. Thus, the length of the vortex trail in- 
creases rapidly and the motion of the vortices becomes dominated by the 
mutual interaction of the vortices rather than by the horseshoe vortices 
at their tail end. The motion of each new vortex ring differs from the 
one shed previously partly because of the continued aging or diffusion 
of vorticity (i.e., a weaker ring is formed at succeeding steps) and 
partly because the motion of each new ring is controlled by all the rings 
shed previously and by the ever-increasing vortex trail. These phenomena 
lead to the rapid demise of vortex rings. In fact, shortly after the 
start of the motion of the lifting surface the horseshoe vortex and the 
vortex rings play a very minor role on the evolution of the trailing 
vortex. 

The trailing vortices are subjected to three known decay mechanisms; 

(i) Aging of vortices due to eddy viscosity; (ii) sinusoidal instability 
leading to the linking of vortices (known as the Crow instability); and 
(iii) vortex breakdown. 

These mechanisms can act independently or in combination. The viscous 
and turbulent diffusion of vorticity lead to the gradual dissipation of 
the trailing vortices. Clearly, this phenomenon is strongly influenced 
by the initial vorticity distribution and the ambient turbulence. Further- 
more, the proximity of other vortices, free surface, rigid surfaces 
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(e.g., the submerged body, ground, etc.), ambient shear, stratification 
of the fluid medium, etc., affect the demise of the vortices. Clearly, 
the aging of the vortices is ever present and affects the evolution of the 
trailing vortices regardless of whether the Crow instability and/or 
vortex bursting take place or not. 

Breaking up by linking or sinusoidal instability involves the growth 
of a rather regular sinusoidal pattern, usually symmetrical for the trail- 
ing vortex pair. This is commonly referred to as Crow linking. The 
time to link (the formation of a series of vortex rings and Horseshoe 
vortices) has been found to depend, statistically, on atmospheric turbu- 
lence within the inertial subrange of eddy sizes (a size range from a few 
cm to some hundreds of meters). The scale of linking varies around nine 
times the vortex spacing. Thus, it is the eddy sizes in the range of one 
to ten times the span or vortex spacing which has the main effect on link- 
ing. The formation of vortex rings by linking (distinctly different from 
those generated at the start of motion) lead to the rapid demise of the 
vortex trail. Consequently, one would want to accelerate the formation 
of links (e.g., by periodically changing the angle of attack of the control 
surfaces) if one is interested in decreasing the hazard posed by the trail- 
ing vortices to other aircraft or in decreasing the energy transferred 
to the fluid medium. 

The vortex breakdown is essentially a single- vortex phenomenon, in the 
sense that it is not necessarily connected with the induced flows asso- 
ciated with the other vortex of the pair. There are, however, frequent 
occurrences of core bursting in conjunction with the sinuous instability, 
in which case the sinuous vortex deformation appears to induce core burst- 
ing. The effect of vortex breakdown is to flatten the velocity profile by 
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significantly reducing the peak velocity (or the velocity of numerous 
peaks of several smaller vortices) under the center of the vortex, but 
without major effect on the speeds away from the vortex center. A 
definite vortex flow remains in all cases, indicating that core bursting 
does not totally eliminate the vortex. In fact, the persistence of the 
outer flow circulation suggests that vortex breakdown may not signifi- 
cantly alleviate the vortex hazard and the energy transfer to the 
surrounding medium. The understanding of the vortex-breakdown and the 
determination of how much circulation remains after breakdown remain as 
critical and vexing questions. 

The trailing vortices move downward (downwash) when generated by a 
lifting surface with a positive angle of attack, as in the case of an 
airplane, or upward (upwash) when generated by a sail plane with a nega- 
tive angle of attack, as in the case of a submarine. In either case, the 
motion of the vortices is determined by the mutual induction, vorticity 
distribution, the three types of the instabilities or demise mechanisms 
cited above, stratification effects of the fluid medium, ambient flow 
turbulence, shear and currents present in the fluid medium, and the 
proximity effects (e.g., free surface, bottom, ground, other bodies, etc.). 
An inviscid flow analysis based on concentrated line vortices can help 
to predict the early stages of motion of the vortices. However, the 
later stages of the motion become increasingly dominated by the factors 
cited above and it is not possible to examine the evolution of the trail- 
ing vortices without taking into consideration the effect of viscosity, 
turbulence, stratification, mutual induction, proximity, and all the other 
phenomena which arise as a consequence of these fundamental variables. 
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It must also be pointed out that whereas the trailing vortices generated 
by an aircraft are under each other's influence (in addition to the 
external non-vortical influences), the trailing vortices generated by a 
large submerged body may also be under the influence of a complex system 
of image vortices imbedded in the body. 

Density stratification within the fluid medium is of great importance 
in the migration and ultimate demise of trailing vortices. The atmos- 
phere is less dense at higher altitude while the oceans exhibit various 
patterns of stratification due to local temperature, salinity, and depth. 
Generally, a vortex migrates with the fluid in which it was born, retain- 
ing most of the original fluid in a so-called recirculation cell. The 
upwash of the vortices in the ocean environment into a layer of lesser 
density gives rise to a buoyancy which is in the opposite direction to 
the velocity resulting from the mutual induction. Furthermore, the motion 
of the recirculation cell produces opposi tely-sigrted vorticity in the 
surrounding medium. This, in turn, gives rise to the entrainment and 
detrainment of fluid from the recirculation cell. Altogether, these 
additional consequences of stratification help to reduce not only the life 
of vortices but also the height to which the vortices rise asymptotically. 
One may also surmise that the stratification will tend to change the 
occurrence and the resulting characteristics of eddy diffusion. Crow 
instability, and core bursting. It is evident from the foregoing that 
the phenomena associated with the generation, migration, interaction, 
and demise of trailing vortices are very complex. As the review of the 
previous investigations will reveal, it is not yet possible to devise 
mathematical models with which some or all of the characteristics of the 



16 



phenomena can be predicted. As in many other flow situations, the role 
played by turbulence and the quantification of the turbulence stresses 
pose difficult questions. The gradual understanding of the behavior of 
trailing vortices must rely for the time being on carefully conducted 
experiments without losing sight of the fact that one must ultimately 
devise a model to acquire a power of prediction as well as a means to 
reduce the number of experiments for model verification. 

It is clear from the foregoing that the reason for developing as much 
an understanding as possible of the behavior of trailing vortices is to 
permit design and improvement of an operational system which properly 
balances safety and efficiency in the motion of vortex-generating bodies 
and also to aid in the development of devices to enhance vortex decay, 
limit vortex upwash, and minimize energy transfer to the surrounding 
medium. The effort towards the understanding of some of the underlying 
phenomena forms the basis of the present investigation. 
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II. REVIEW OF PREVIOUS INVESTIGATIONS 



Vortices and vortex wakes have become a major theme of aerodynamics 
research since the advent of the large aircraft. The objectives have been 
to predict the trajectories and lifetimes of trailing vortices, in order 
to relieve the constraint on flight operations imposed by the vortex 
hazard. Much of the progress made during the past two decades was dis- 
cussed at the Symposium on Aircraft Wake Turbulence and Its Detection 
[Ref. 1], held in Seattle in 1970, and at the Aircraft Wake Vortices 
Conference [Ref. 2], held in Cambridge, MA in 1977. 

The fundamental questions regarding the evolution of the trailing 
vortices may be summarized as follows: (i) how does the vorticity gener- 
ated by the lifting surfaces organize itself into two or more centers of 
vorticity behind and to each side of the lifting surface? (ii) how do 
the discrete centers of vorticity interact and affect the stability of 
each other? (iii) how does a vortex demise and how does the demise of a 
single vortex differ from that of a vortex pair? And finally, (iv) what 
is the effect of stable stratification on the generation, roll-up, inter- 
action, stability, and demise of vortices? 

A. ROLL-UP OF VORTICITY 

It is a well known fact that a vortex sheet is trailed from a finite- 
aspect-ratio lifting surface as a consequence of the nonuniform spanwise 
lift on the wing. The roll-up of this sheet into discrete vortices as 
a result of a convective motion was recognized as early as 1907 by 
Lanchester [Ref. 3]. Since then there have been numerous studies to 
describe this complicated phenomenon. The first mathematical model to 
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describe the inviscid structure of the vortex wake was proposed by Prandtl 
[Ref. 4], He based his model on the conservation of mechanical energy 
and found that the vortex tube radius for an elliptically loaded wing 
is approximately 8 percent of the wing span. In this analysis the choice 
of the swirl velocity distribution is somewhat arbitrary, being only 
constrained by circulation at large radius from the vortex center and 
the integral constraint on kinetic energy. In other words, the model 
does not address itself to the determination of the final structure of 
the wake for a given initial vorticity distribution. 

The second model, which overcomes some of the arbitrariness inherent 
in the Prandtl model, was developed by Betz [Ref. 5]. Betz's model is 
based on four integral invariants of an incompressible two-dimensional 
bounded vorticity distribution. The Betz model has now been checked 
against measurements made behind model and full-scale aircraft and it 
has been found that the model predicts distributions which are irf far 
better agreement than the Prandtl model [Ref. 4]. 

In more recent times the advent of numerical methods and computer 
gave impetus to the use of discrete vortex models. In most numerical 
analyses, the attached boundary layers on the wing surfaces are repre- 
sented by bound- vortex sheets while the separated free-shear layers are 
modeled by free-vortex sheets. The free- vortex sheets join the bound- 
vortex sheets along the separation lines which are known a priori for 
wings with sharp edges. Furthermore, it is assumed that the vortex 
breakdown (see e.g.. Refs. 6-8) occurs far downstream so that the 
primary core size and its variation in the vicinity of the wing are 
negligible. This assumption limits the angle of attack (dependent on 
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the aspect ratio of the wing) at which the inviscid model is applicable. 
Furthermore, the flow outside the bound-vortex sheet(s) and the free- 
vortex sheet(s) is assumed irrotational , Within these limitations, the 
resulting potential flow model represents the main features of the real 
flow to a high degree of accuracy. The literature contains several steady 
and unsteady inviscid-flow models with various degrees of limitations 
and drawbacks. 

The first group of models uses slender-body theory and conical flow 
assumption [Refs. 9-12]. These models satisfactorily predict the pressure 
distribution over the front portion of the wing surface. In the rear 
portion, the models fail to predict a satisfactory pressure distribution 
because the Kutta-Joukovsky condition is violated at the trailing edge. 
Such models are limited to simple delta planforms. The second group of 
models uses a nonlinear discrete-vortex method [Refs. 13-20]. In these 
models, the bound vortex sheet and the free-vortex sheets are approximated 
by a set of concentrated vortex lines. The bound vortex sheet is replaced 
by a bound-vortex lattice, while the free-vortex sheet is replaced by 
segmented free-vortex lines (in the case of steady flow) or by growing 
free-vortex lattice (in the case of unsteady flow). The boundary condi- 
tions are satisfied at certain control points on the bound- and free- 
vortex system using an iterative technique. 

Although the discrete vortex model has worked well for many years for 
attached flow problems, when vortex-type separation from leading edges 
and/or tips occurs, the free-vortex system lies close to the lifting- 
surface bound vortices and the results are found to be sensitive to 
variations in the shapes of the vortex elements and the relative lengths 
of the vortex segments (see e.g.. Ref. 21). 
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The third group of models employs doublet-panels (see e.g.. Refs. 
22-23). In this method, the wing and its free-vortex sheets are divided 
into networks of quadrilateral panels. Each panel of the networks 
representing the wing has a biquadratic local doublet distribution and 
a bilinear local source distribution. Source and doublet splines are 
used to express the distributions of singularities on the networks in terms 
of discrete values of singularity strength at certain standard points on 
each network. The boundary conditions and continuity of singularity 
strengths across abutting networks are enforced at certain standard points 
on each network. The results of this method are generally good when the 
solution converges (see e.g.. Refs. 24-25). Apparently, the difficulty 
in obtaining convergence is due to the failure in satisfying the con- 
tinuity of the derivatives of the doublet strength across abutting net- 
works. This is equivalent to the existence of concentrated vortex lines 
between abutting networks. Additional related works are discussed in 
Ref. 26 in connection with the AGARD Symposium on High Angle of Attack 
Aerodynamics . 

The fourth group of models employs a nonlinear hybrid vortex method 
(see e.g.. Refs. 27-28). In this method, continuous-vorticity and vortex- 
line representations of the wing and its separated free-shear layers are 
used. Continuous vorticity is used in the near-field calculations while 
discrete vortex lines are used in the far-field calculations. The calcu- 
lated net surface pressure and section normal-force coefficients for 
wings with tip separation are in good agreement with the experimental 
data and the nonlinear discrete-vortex method. 

The roll-up calculations do not give details of the structure of the 
tip vortex to be used as initial conditions in a viscous solution. Often 
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the approximate solution of Kaden [Ref. 29] for the unsteady roll-up of 
an infinite sheet is used as the basis for the prediction of the decay 
of laminar trailing vortices [Ref. 30]. The most obvious drawback of 
the existing models is the lack of a realistic model of the leading- 
edge core and its feeding sheet. Modeling of the leading-edge core and 
its feeding sheet had been first introduced by Brown and Michael [Ref. 
31]. This model has been improved by Mangier and Smith [Ref. 32] and 
Smith [Ref. 33]. Additional improvements have been made by Moore [Ref. 
34] who devised a method to collapse the tightly-wound inner portion of 
the vortex sheet into a single tip vortex. This was done primarily to 
avoid the numerical difficulties associated with the crossing of the 
vortex sheets and with the spurious interaction between the successive 
turns of the inner portions of the vortex sheet. Clearly, this sacri- 
fices the detail of the tip vortex structure. The latest efforts regard- 
ing the calculation of the roll-up process use the method of rediscreti- 
zation of shear layers. This method has been devised by Fink and Soh 
[Ref. 35] and extensively used by Sarpkaya and Shoaff [Ref. 36]. The 
numerical models based on the inviscid flow theory do not account for 
the diffusive effects of viscosity and turbulence. Thus, they are not 
capable of predicting the evolution of the trailing vortices. 

B. TRAILING VORTICES 

The vortex sheets roll-up asymptotically and it may be several wing 
spans behind the aircraft before the fully rolled-up condition is 
reached. With the high Reynolds numbers encountered on lifting surfaces 
in free flight, it is likely that their boundary layers, and thus their 
trailing vortex sheets, are turbulent. At finite Reynolds numbers. 
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vorticity diffuses, and the sheet becomes a vortex layer of finite, non- 
zero thickness. Successive turns in the centermost portions of the 
rolled-up layer merge, the spiral structure disappears and, if the radius 
of the vortex core is greater than the diffusive length scale, an equi- 
librium structure ensues. Thus, one anticipates a multi-structured 
core, the innermost part of which has effectively reached a state of 
equilibrium and surrounding this an annulus in which the turns have merged 
or are in various stages of merging, and an outer region, beyond which 
discrete turns of the spiral remain evident. 

Dynamics of the flow require axial velocities to be present in trailing 
vortices [Ref. 37] and these can have an important effect on the stability 
of the trailing vortex systems and play an important role in the occurrence 
of vortex breakdown. It is also a well-known fact that the occurrence of 
vortex breakdown in trailing vortices may be delayed or enhanced by the 
acceleration or deceleration of the ambient flow due to the proximity of 
streamlined axisymmetric bodies whose control surfaces may have given rise 
to the development of trailing vortices. Thus, it is important to make 
a distinction between the trailing vortices of an aircraft and those of a 
submerged body. In the case of the former, the trailing vortices are not 
affected by the proximity of the aircraft body partly because the length 
of the aircraft in terms of the wing span is relatively small and partly 
because the wing span is many times that of the body diameter. In the 
case of a common submerged body, the wing span of the major control 
surfaces is in the order of the body diameter and the length of the body 
in terms of the wing span may be in the order of ten. Consequently, the 
trailing vortices of a submerged body are significantly affected by the 
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flow about the body and by the propulsion system which may modify the flow 
about the body so as to affect the evolution of the trailing vortices 
and the occurrence of vortex breakdown. Evidently, the said differences 
stem from the fact that the wings on an aircraft provide the lift neces- 
sary for flight whereas the control surfaces on a submerged body provide 
only partial lift for depth and attitude control. The major portion of 
the lift is already provided by the buoyancy of the body. 

Numerous attempts have been made to model the behavior of laminar and 
turbulent trailing vortices. Direct calculation of the structure of the 
trailing vortices is done by assuming a final flow configuration and re- 
quiring conservation of some chosen flow quantities. Most early discus- 
sions of the far-field structure were based on the vortex -wake model of 
Spreiter and Sacks [Ref. 38]. They have assumed that the wake consists 
of two Rankine vortices, with uniform vorticity inside the core and zero 
vortici1:y outside. The circulation and vortex separation are given by 
Stokes' theorem and conservation of momentum. The core radius is deter- 
mined by requiring kinetic energy to be conserved. For an elliptically 
loaded wing this model gives a core radius of about 0.1 times the trail- 
ing vortex separation. Test flight data reported that the cores are 
somewhat smaller than predicted and that all of the vorticity is not in 
the core [Refs. 1-2]. 

Donaldson [Ref. 39] noted that the model of Betz [Ref. 5] which assumes 
that each half of the vortex sheet rolls up into a circular vortex core 
conserving circulation, impulse, and momentum of impulse, could be used. 
This conserves the moment of inertia of the vortex distribution about its 
centroid. Betz further determined the distribution of circulation in the 
vortex by assuming the vorticity outboard of some point, which will form 
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the core of the vortex inside some radius, conserves its moment of iner- 
tia about its own centroid during roll-up. Both Jordan [Ref. 40] and 
Donaldson [Ref. 41] have obtained similar solutions. This approach 
agrees much better with flight test data than the model of Spreiter and 
Sacks [Ref. 38]. 

C. STABILITY OF TRAILING VORTICES AND THE EFFECT OF STRATIFICATION 

If the trailing vortices could be represented as a two-dimensional 
vortex pair, the study of their motion may be more straightforward. 
Unfortunately, they are most unstable to three-dimensional disturbances 
so that the vortices rapidly lose their two-dimensional character. This 
is one of the major differences between the motion of a single vortex 
and the vortex pair. A single vortex is known to be very stable whereas 
a vortex pair suffers from numerous instabilities. Nevertheless, the 
study of two-dimensional vortex pairs has continued in hopes of obtaining 
a satisfactory solution that will match what little three-dimensional 
data are available. The majority of the data have been obtained in large 
scale atmospheric tests using aircraft [Refs. 42-43] or from vortex rings 
[Refs. 44-45]. Barker and Crow [Ref. 46] and Tomassian [Ref. 47] have 
created two-dimensional vortices in a small tank with large circulations 
and comparably small separations. There are significant differences 
between the experimental results cited. These differences and their 
explanations are further compounded by the contradictions of the analyses. 

There are fundamental differences in the assumptions made and widely 
differing conclusions have been reached regarding the motion of the trail- 
ing vortex pair in a stratified medium. The lack of consistent laboratory 
data did not help to resolve the differences and to delineate the most 
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important elements of the physics of the motion. Most models use an oval 
shaped fluid body carried by the vortices under the influence of their 
mutual induction and the density difference with the surroundings. There 
is a significant difference in the accounting of fluid momentum, impulse, 
energy, viscous and turbulent diffusion, and in the assumptions concern- 
ing the buoyancy induced vorticity entrainment and detrainment within 
the recirculation cell (see Fig. 1). One school of thought holds that the 
recirculation cell buoys upward in a stably Stratified atmosphere [Refs. 
48-49], while others maintain that the vortices accelerate downward 
[Refs. 50-51]. The resolution of the dilemma is important, because of 
the understanding of wake hazard and of the energy transferred to the 
surrounding medium. 

Table I, adapted from Widnall [Ref. 52], describes the differences in 
some of these models. The models of Scorer and Davenport [Ref. 50] and 
Crow [Ref. 51] do not include turbulence effects. The models of Hecht 
et al . , [Refs. 44 and 53] are computer generated flow fields assuming 
constant recirculation cell areas. Clearly, the differences in these 
models cannot be resolved until reliable laboratory data are available. 

It is important that the experiments model the trailing vortices and not 
just a pair of two-dimensional vortices. The three-dimensional nature 
of the flow and the importance of the axial velocity in the trailing 
vortices must be reflected in the proper simulation of the behavior of 
vortices. It is, thus, important that the differences between the three- 
dimensional experiments and the two-dimensional analytical or numerical 
models be clearly delineated. It is only through such comparisons that 
one would be able to determine as to whether enough physics have been 
incorporated into the analytical model. 
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Figure 1. Downward Migrating Recirculation Cell 
in a Stratified Atmosphere. 
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TABLE I 



Comparison of Theoretical Studies of Vortex 
Pairs in a Stably Stratified Atmosphere 



Author 


Vortex 

Spacing 


Migration 

Speed 


Buoyancy 

Generated 

Vorticity 


Costen [52] 








Weak stability 


Decreases 


Increases 


No Effect 


Strong stability 


Increases 


Stops 


- 


Kuhn & Nielsen [58] 


Decreases 


Increases 


Partly Ent. 


Saffman [49] 


Increases 


Stops 


Entrained 


Scorer & Davenport [50] 


Decreases 


Increases 


Detrained 


Tombach [59] 








Strong stability 


Increases 


Stops 


Entrai ned 


Weak stability 


Decreases 


Increases 


Entrained 


Tulin & Shwartz [48] 


Increases 


Stops 


Entrained & 
Annihilated 


Crow [51] 


Decreases 


Increases 


Created in 
Wake 


Narain & Uberoi [60] 


Decreases 


Decreases 


Entrai ned 


Hecht et al . [44] 


Decreases 


Stops 


Created at 
Boundary 


Hecht, Bilanin, 
Hirsh [53] 


Not 

Appreciable 


Stops 


Created at 
Boundary 
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The first quantitative analysis of the three-dimensional sinusoidal 
instability observed in trailing vortices in an ideal homogeneous fluid 
was given by Crow [Ref. 54], His results showed good agreement with the 
observations of the general features of aircraft vortex instability. 

Crow and Bate [Ref. 55] explored the growth of vortex-pair instability 
when driven by ambient turbulence. They have calculated the time to 
linkage as a function of the Kolmogorov turbulence parameter. The possi- 
bility of exciting the vortex-pair instability by changing the lifting 
surface span-load distribution was suggested by Crow and demonstrated in 
flight by Chevalier [Ref. 56] by changing aircraft pitch angle periodically. 

The trailing vortices may also suffer from instabilities inherent 
to single vortices, i.e,, from vortex breakdown. Although the mechanism 
is not well understood, the presence of axial velocities in the trailing 
vortices and the proximity of the body from which the vortices are 
emanating play significant roles in the occurrence of vortex breakdown. 

The vortex breakdown has been studied in detail for single vortices 
(see e.g., Refs. 6-8). It is not known what would be the effect of mutual 
induction and straining of vortices on the breaking of one or both of 
the vortices. Furthermore, it is not clear as to what happens after the 
breakdown. It was assumed that the breakdown leads to the total dissi- 
pation of circulation and thus its occurrence is most desirable. However, 
observations [Ref. 2] have shown that the vortex breakdown does not 
totally dissipate the existing circulation. Thus, the determination of 
how much circulation remains in a vortex after the breakdown constitutes 
one of the major questions. One would also like to determine as to how 
long the remaining circulation persists. Qualitatively, it is now clear 
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that the path to the ultimate destruction of trailing vortices has not 
yet been charted. The Crow instability, laminar and turbulent diffusion, 
entrainment and detrainment of oppositely-signed vorticity, vortex break- 
down or a cascade of vortex breakdowns, and the proximity effects (body 
or ground) play significant roles. For a submerged body, the deceleration 
of the ambient flow toward the tail of the body may give rise to vortex 
breakdown at a fixed point relative to the body. This, in turn, results 
in a series of vortex breakdowns as the body moves at a constant speed, 
leaving behind a cascade of breakdowns. Since each breakdown separates 
a supercritical swirling flow from a subcritical flow (downstream of the 
breakdown) (see Ref. 7) and since the subcritical flow cannot sustain 
waves, each breakdown rapidly moves upstream to catch up with the newest 
breakdown. This, in turn, leads to rapid weakening of the trailing vortices. 
The quantification of these qualitative observations constitute the essence 
of future research efforts. 

D. SUMMARY OF THE PREVIOUS INVESTIGATIONS 

It is clear from the foregoing that a great deal of additional work 
is needed to understand the rolling-up of vortices in both the strati- 
fied and unstratified fluid medium, the evolution of trailing vortices, 
the demise mechanisms, and the processes leading to the ultimate destruc- 
tion of vortices. Analytical studies must be based on reliable experimental 
data in order to incorporate as much relevant physics as possible into 
the models. Furthermore, the analyses must be based on the observations 
and measurements of three-dimensional trailing vortices rather than those 
of two-dimensional vortices in order to bring into focus the important 
role played by the three-dimensionality of the flow and the axial velocity 
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in the vortices. In fact, a rigorous investigation of the entire phenomenon 
might begin with trailing vortices generated by relatively simple bodies 
(such as a Delta wing or rectangular wing) in an unstratified medium and 
then expanded to include the effect of stratification. At present there 
are no extensive experiments in either medium and no detailed calcula- 
tions have been reported. 
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Ill, EXPERIMENTAL EQUIPMENT AND PROCEDURES 



It is clear from the foregoing that there is very little or no data 
regarding the motion of trailing vortices in stratified and unstratified 
fluids. It is also clear that the analytical models and numerical calcu- 
lations can be judged, validated, and upgraded only through the use of 
data obtained with three-dimensional trailing vortices under controlled 
laboratory conditions. This, in turn, requires the use of a large water 
basin in which various models may be towed at desired speeds. 

The test facility used in the present investigation was developed 
with the following specific objectives in mind: (i) the establishment of 
the parameters governing the phenomenon and the scaling laws; (ii) the 
understanding of the characteristics of the sinusoidal instability in three- 
dimensional trailing vortices with axial velocity; (iii) the determination 
of the role played by the vortex breakdown or cascade of breakdowns and of 
the circulation remaining in the trailing vortices following breakdown; 

(iv) the establishment of any interaction between the sinusoidal instability 
and the vortex breakdown; (v) the investigation of the effect of ambient 
stratification, ambient turbulence, and model-induced turbulence on the 
evolution and demise of trailing vortices; (vi) the determination of the 
characteristics of internal waves generated by the body and/or trailing 
vortices prior to, during, and after the onset of sinusoidal instability 
and/or vortex breakdown; (vii) to examine the behavior of trailing vortices 
generated by a lifting surface undergoing time-dependent oscillations (with 
or without the body- proximity effects); (viii) to study the evolution of 
Horseshoe vortices and vortex rings; (ix) to examine the effect of the 
alterations in the axial velocity of the vortices at the time of their 
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generation; and (x) to investigate the interaction of the trailing vortices 
with the vortices generated by the body and/or control surfaces. 

The experimental equipment consists of the water basin, the strati- 
fication system, and the test models. 

Ao WATER BASIN 

The basin is 36 feet long, 3 feet wide and 5 feet deep and is made of 
0.5-inch thick aluminum plates, welded appropriately (see Fig, 2), The 
back wall of the tank is reinforced to prevent bulging under hydrostatic 
pressure. The front wall is fitted with square windows of 1,0-inch thick 
plexiglass for flow visualization. 

Vertical and horizontal scales of 0.5-inch increments are marked on 
the central window. The scales are used during the experiments to evaluate 
the position of the trailing vortices (see Fig. 3). 

The aluminum interior of the tank is painted with marine epoxy to 
prevent corrosion. In addition, the inside face of the back wall is painted 
with flat white paint to provide better contrast and make the flow visuali- 
zation easier with the dye used. 

The filling pipes lie on the bottom of the tank. They are plastic 
pipes 2,0-inch diameter, perforated along their entire length with holes 
of 1/16-inch diameter (see Fig. 4). 

Two drains are provided at the bottom at each end of the basin. Two 
parallel rails of 1.0-inch diameter are mounted along the bottom of the 
tank, A carriage rides smoothly on these rails and provides the test 
body with a constant velocity through the use of an endless cable and a 
3/4 Hp DC motor (see Fig, 5). The velocity of the model is measured and 
continuously monitored through the use of a magnetic velocity sensor. 
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Figure 2. Water Basin. 
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Figure 3. Central Window and the Horizontal and Vertical Scales. 





36 



Figure 4. Top View of Stratification System. 
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Fi'jure 5, Towing Mechanism. 



mounted above a precision gear. The sensor counts the number of teeth 
and yields and instantaneous velocity proportional to the velocity of 
the model . 

The two rails, the carriage and the filling pipes are located on or 
near the bottom of the basin and under a turbulence management system. 

This system consists of a one-inch-thick polyurethane foam, sandwiched 
between two perforated aluminum plates (see Fig. 6). 

B. STRATIFICATION SYSTEM 

The stratification system consists of four tanks (60-gallon capacity 
each) and the necessary tubes and valves. 

First, brine water of certain salt concentration is prepared in one 
of these tanks. Then precalculated amounts of brine water are added to 
the other three tanks, depending on the desired density gradient and the 
number of the layers. The calculation of the amount of salt needed, for 
a certain density gradient or for a certain stratification parameter, is 
discussed in Appendix A. 

The resulting solution is drained into the basin at a rate of 2.1 CFM. 
Necessary precautions have been taken to prevent air entrainment and 
subsequent bubble generation. Furthermore, two air-bubble traps have 
been, installed, one just before the main valve of the filling system 
and one at the middle point of the delivery pipes. 

This system can fill the water basin and establish a desired stratifi- 
cation in about three hours. There is no difficulty regarding the desired 

C 

density gradient, even for the lower gradients (e.g., dp/dy = 5x10’ 
slug/ft^/inch) . The density discontinuity between the 0.9- inch layers 
disappeared by diffusion after a short time for the stratifications of 
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Figure 6. Turbulence Management System. 



interest. Smooth density gradients were obtained as verified by a 
salinity meter. The parts of the stratification system are shown in 
Figs. 4 and 7. 

C. MODELS 

Three models were used in the experiments. The first two models are 
Delta wings, made of 1/8-inch thick aluminum with hollow interiors. 

Each model cavity is filled with dye prior to the filling of the basin 
with water. The two holes at the rear edge of each model and the third 
hole in the middle of the upper surface are plugged prior to an experi- 
mental run to prevent dye leakage (see Fig. 8). The dye consisted of 
food coloring and proper amounts of alcohol. The density of the mixture 
was such that the dye was only extremely slightly buoyant relative to 
its surroundings . 

The models are mounted on their bases by means of a streamlined thin 
aluminar bar (its cross section was a NACA 0006 foil) and set at the de- 
sired angle of attack. The models are placed at a level 9 inches above 
the foam. The first Delta wing had a vertex angle of 39.5055 degrees and 
a base length of 7.9 inches. The second model had a vertex angle of 50 
degrees and a base length of 5.071 inches. 

The third model is a submerged planing surface made of plexiglass and 
attached to a streamlined body. The body is cut in half longitudinally 
with the bottom half removed to allow testing at greater depths with 
negligible hydrodynamic ground effects. The planing surface has a base 
length of 6.96 inches with an aspect ratio of 2.20 (Aspect ratio = [base 
length] /Area). The model is mounted to the carriage such that the plane 
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Figure 7. Stratification System Auxiliary Tanks. 






Delta-wing Model No. 1 



DIMENSIONS 



Del ta 1 : 



B = 7.90" 

b = 5.585" 

0 

h = 11 .00" 

0 

e = 39.5 deg 



Delta 2 : 

B = 5.07" 

b = 3.70" 

0 

h = 6.00" 

0 

e = 50.0 deg 




Figure 8. Delta-Wing Models. 
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bottom of the body is flush with the top of the aluminum plate above the 
foam. 

A hollowed core within the supporting column of the planing surfaces 
is filled with dye. The dye is dispensed from the tips of the planing 
surfaces with small holes connected to the dye reservoir (i.e., the 
hollow core) . 

As noted earlier all models are pulled by means of a DC motor, pulley, 
and cable system at the desired speed (ranging from 1.0 to 5.0 ft/sec). 

The entire arrangement is shown in Figs. 9 and 10. 

D. TEST PROCEDURES 

The model is initially filled with dye and all the fill and drain 
holes on the model are closed with wooden pegs. Then the basin is filled 
gradually either with fresh water (for unstratified flow experiments) or 
with salt water (for stratified flow experiments) to the desired level. 
After sufficient period of waiting for the elimination of any internal 
currents in the basin the model was set in motion at a desired speed. 

The motion of the trailing vortices are recorded on high-speed Polaroid 
film at the test section (one of the plexiglass panels near the middle 
of the basin). Each picture included two clocks accurate to 0.01 second, 
the vertical and horizontal scales on the plexiglass window and, of 
course, the side view of the trailing vortices as they rose from the model 
after formation (see Fig. 11). The time interval between successive pic- 
tures is determined from the two clocks. The vertical rise of the vortices 
is determined from the vertical scale. Attention has been paid to the 
fact that the vortices are farther away from the scale on the plexiglass 
window and that the scale placed vertically in the middle of the test 
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Figure 9. Motion Generation Mechanism. 
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Figure 10. Front View of the Basin. 
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Figure 11. Side View of the Trailing Vortices. 





section does not exactly correspond to the scale marked on the plexi- 
glass window due to refraction and parallax. The necessary correction 
was made by photographing a scale placed in water in the middle of the 
test section together with the scale marked on the plexiglass window. 

This resulted in a simple conversion table which enabled the determina- 
tion of the actual position of the vortex from the scale reading on the 
photograph . 

The results are non-dimensional ized and plotted in various forms and 
compared with those obtained in the previous runs. Each experiment was 
repeated at least twice to ascertain that the experiments were repeatable. 

The experiments with the formation of Horseshoe vortices and vortex 
rings were conducted in the manner described above with the exception that 
the model was initially placed closer to the test section, rather than 
at one end of the basin, and that both the horizontal and vertical scales 
were used to determine the velocity components of the translation of the 
vortices . 

All trailing vortices and vortex rings were recorded on film until the 
time they have completely dissipated either due to aging (diffusion of 
vorticity due to viscosity, turbulence, entrainment, and detrainment) , or 
due to Crow instability (sinusoidal instability and linking leading to 
the formation of vortex rings), and/or due to vortex breakdown (core 
bursting). It was thus possible to determine the life span of vortex 
cores from their formation to their final demise. 
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IV. RESULTS 



A. DIMENSIONAL ANALYSIS 

The experimental procedure is based on a careful dimensional analy- 
sis of the parameters entering into the phenomenon. 

The dependent parameter of major importance is the displacement H 
of the vortex in the vertical direction and it may be expressed as a 
function of the following independent parameters for a given body shape: 



where t denotes the time; U, the forward velocity of the model; D, the 
initial depth of submergence of the vortex pair; p^, the fluid density; 
9p/3y, the density gradient; v, the kinematic viscosity of water; B, the 



base width of the model (see Fig. 8); e, the vertex angle; a, the angle of 
attack; and g, the gravitational acceleration. The other parameters such 
as the vortex strength r, the core radius of the vortices r^, the vortex 
separation b^, and the velocity of vertical migration of vortices V^, 
are determined by the same parameters governing the variation of H. 

A dimensional analysis of Eq. (1) yields 




( 1 ) 




( 2 ) 



in which 




1/2 
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and is known as the Brunt-Vaisal a frequency. The other dimensionless 
parameters appearing in Eq. (2) are self-explanatory. 

In the present investigation B, e, a, and 0 were determined by the 
model being investigated. The velocity of the model was varied from 
0.32 ft/s to 4.04 ft/s in experiments with fresh water (N =0). The 
experiments with stratified water were carried out by varying N 
(N = 0.0974 s ^ to 0.2101 s'^ ) at various velocities. 

The dimensionless parameter NB/U is often replaced by another parameter 
Nb^/V^, known as the stratification parameter. It should be noted that 
both b^ and are dependent parameters, b^ is determined by the shape of 
the body and its angle of attack. is determined by 



h = fr 2 ^ u£ UB X 
U ’^^B’ U ’gB’ V 



(3) 



Thus, Eq. (2) may be recast as 




y t ^ Hb T 

^ ’b^’ ’gb„^ 

0 0 0 ^0 



V b 
0 0 






(4) 



Of the parameters appearing in Eq. (4) only t, D, N, g, and v are inde- 
pendent and may be varied irrespective of the shape of the body. However, 
the remaining parameters with the exception of a depend on the shape, 
aspect ratio, and the velocity of the body. The variation of e and a 
are imbedded in the variation of and b^. Experiments were performed 
with a given body by varying the angle of attack and the speed of the 
body for a given Brunt-Vaisal a frequency, N. Then and H were deter- 
mined from the pictures. The parameters s and b^ were obtained from the 
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existing experimental data [Ref. 61]. Subsequently, the data were plotted 
in terms of H/b^ versus V^t/b^ for various values of D/b^, and 

2 

The measurements have shown that the dimensionless parameters V^/gb^ and 
are not important within the range of the parameters encountered 
in the experiments. With the foregoing arguments, Eq. (4) may be written 
as 




f( 





(5) 



This equation formed the basis for the analysis of the entire experimental 
data . 



B. DELTA-WING GENERATED VORTICES IN UNSTRATIFIED WATER 

Sample data obtained with the first Delta-wing (referred to hereafter 
as Delta 1) in unstratified water (Nb^/V^ = SP = 0) are shown in Fig. 12 
for two values of and for D/b^ = 5.19. Two symbols, identifying two 
independent runs, show that there is very little scatter in the data. Had 
the vortices maintained their strength and spacing, as they would in an 
ideal fluid, H/b^ would have been equal to V^t/b^. This straight line 
relationship is also shown in Fig. 12. Clearly, all the data points are 
below the ideal line and show that the vortices do not retain their original 
characteristics. One may argue that this may be as a consequence of 
several factors such as the decrease of circulation, the increase of the 
vortex spacing b, and the proximity of the free surface. The variation of 
b with time may be excluded from this discussion because the observations 
and measurements have shown that it remains practically constant for large 
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enough values of D/b^. As to the free surface effects, it will be shown 
that they do not play any role in the rise of the vortices for D/b^ 
larger than about 4.3. Thus one may conclude that for large D/b^ values, 
the difference between the actual and the ideal values of H/b^ is due to 
the decrease of the circulation of vortices with time. This decrease, as 
noted in connection with the discussion of previous studies, is a conse- 
quence of eddy diffusion of vorticity. Crow instability, and vortex break- 
down. This fact has not been recognized in the past because of lack of 
reliable data and it was often assumed that the vortex strength remained 
constant. This in turn gave rise to conflicting conclusions cited in 
Table I. 

Figure 13 shows a composite plot of the data obtained with Delta 1 

for various values of V at a fixed value of D/b„ = 4.31. The primary 

purpose of this figure is to substantiate the assertion that the dimension- 
2 

less parameters V^/gb^ (a representati ve Froude number) and V^b^/v (a 

representative Reynolds Number) do not play experimentally detectable 

roles in the variations of H/b with V t/b , at least for unstratified 

0 0 0 

flow experiments. The two-fold variation of for constant values of b^, 
g, and v yields a four- fold variation in the Froude number and a two-fold 
variation in the Reynolds number. Clearly, Fig. 13 does not exhibit such 
a dependence on the said two parameters. In fact, it is because of this 
finding that the governing parameters were finally reduced to that given 
by Eq. (5). 

The next question to be addressed is the dependence of H/b^ on D/b^. 

In other words, at what D/b^ does the free surface affect the rise of the 
vortices and in what form does this effect exhibit itself? 
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Figure 13. Vertical Rise for Delta 



Figure 14 shows the data obtained with D/b^ values ranging from 2.15 
to 5,19, again in unstratified water. Evidently, the variation of H/b^ 
versus V^t/b^ for D/b^ = 4.31 and D/b^ = 5,19 is such that the differ- 
ences are negligible and certainly well within the experimental errors. 

For other values of D/b^, the vortices reach their limiting height very 
quickly and then begin to bounce back from the free surface. The smaller 
the D/b^, the smaller is the limiting height. These experimental facts 
may be explained partly in terms of the behavior of inviscid vortices in 
the vicinity of a plane wall. 

The free surface acts as a reflection plane and the motion of the 
vortex pair is determined not by their mutual induction alone but also by 
the image vortex pair above the free surface. In other words, the tra- 
jectory of the ideal vortices is determined by the mutual induction of 
four vortices (see Fig, 15).. Lamb [Ref. 62] has shown that the trajectories 
of line vortices in ideal unstratified fluid are hyperbolas given by 



(2x/b^)' 



(2y/b^)' 



= 1 + 



b , 
(— )' 



( 6 ) 



where b^ is the initial spacing of the two real vortices, and D is their 

initial depth. Equation (6) shows that the vortices move sideways as 

they approach the free surface. For large values of x and D, y approaches 

half the initial separation, i.e., the vortex pair does not come closer 

to the free surface than a distance of b /2. 

0 

A similar analysis cannot yet be carried out in a real fluid including 
the effects of stratification and turbulence. The strength of the vortices 
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Figure 14. D/b Effects on Vertical Vortex Rise for Delta 




Figure 15. Mutual Induction of Four Vortices and the 

Free Surface Effects. 
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and their rise velocity continue to decrease (see the slope of the data 
for D/b^ = 2,15 in Fig. 14). 

As the vortices move too far sideways, they begin to fall under the 
influence of additional image vortices due to the presence of the side 
walls. Even though it is not a matter of primary importance, one can 
show that the combined effect of the image vortices due to the side walls 
and the free surface is to push the vortices downward. It is not yet 
possible to verify whether there exists a physical mechanism which will 
make the vortices bounce from the free surface had there been no side 
wall effects. It is possible that the decrease in the mutual influence 
of the trailing vortices may give rise to alterations in the velocity 
fields in the presence of a movable free-surface boundary and this in 
turn may result in a bouncing back of the vortices. In any case, only 
the experiments in a large basin can resolve the question. 

An additional fact which can be extracted from Fig. 14 is the lack 
of importance of the Froude and Reynolds numbers on the rise of vortices 
even when the free surface effects predominate. The data for all D/b^ 
values were obtained at different model speeds even though they have been 
shown with a single symbol in Fig. 14. The plotted as well as the tabu- 
lated data (see Appendix C) show that (see for example the data for 
D/b^ = 2.15 in Fig, 14) neither the Froude number nor the Reynolds number 
plays any role on the rise of vortices in unstratified flow with or 
without free surface effects, within the range of the parameters encoun- 
tered in this investigation. 

The results presented above were obtained with Delta 1. In order to 
ascertain further the validity of the conclusions reached regarding the 
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role played by the governing parameters (see Eq. (5)) experiments were 
carried out with the second Delta-wing model (Delta 2). Figure 16 shows 
H/b versus V t/b^ for D/b = 7.73 in unstratified water for various 
values of (Tabulated data are presented in Appendix D) . Clearly, the 
results do not depend on and hence on the model speed. The comparison 
of the data obtained with the two models is presented in Fig. 17. Aside 
from some minor scatter at large values of V^t/b^, the results of the 
two models compare quite well. Thus, it may be concluded on the basis 
of the extensive data collected from the two models that in unstratified 
water the rise of the trailing vortices is governed by the dimensionless 
parameters H/b^, V^t/b^, and D/b^. In all cases the vortices gradually 
lose their strength and thus they move upwards with ever decreasing 
velocities until they are dissipated by various demise mechanisms dis- 
cussed earl ier . 

C. DELTA-WING GENERATED VORTICES IN STRATIFIED WATER 

Experiments in stratified water were carried out with the same Delta- 
wing models for various stratification parameters (SP = Nb^/V^). First, 
the results obtained with Delta 1 for five SP values will be presented. 
Then the results obtained with Delta 2 for three SP values will be 
shown. Then the results obtained with the two models will be compared. 

Figures 18 through 22 show H/b^ versus V^t/b^ for SP values of 0.375, 
0.50, 0.625, 0.75, and 1.0 for a fixed value of D/b^ = 4.31. It is 
important to note that this particular value of D/b^ was chosen so as to 
exclude the effects of free surface and thus to concentrate only on the 
effects of linear stratification. A complete plot of the aforementioned 
data is presented in Fig. 23. Clearly, the larger the stratification 
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Figure 16. Vertical Vortex Rise for Delta 2, SP = 0, D/b = 7.73. 
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Figure 17. Comparison of the Vortex Rise for Delta 1 and Delta 2 for SP 
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Figure 18. Vertical Vortex Rise for Delta 1 , SP = 0.375, D/b 
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Figure 19. Vertical Vortex Rise for Delta 1 , SP = 0.50, D/b 
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Figure 20. Vertical Vortex Rise for Delta 1 , SP = 0.625, D/b = 4.31. 
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Figure 2K Vertical Vortex Rise for Delta 1 , SP = 0^75, D/b = 4.31 
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Figure 22„ Vertical Vortex Rise for Delta 1, SP =U0, D/b = 4.31 
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Figure 23„ Stratification Effects on Vortex Rise for Delta 



parameter, the smaller the rise of vortices. The last data points (of 
each particular run) correspond to the state where the vortex core was no 
longer discernible and the surrounding dye was randomly diffused with no 
sign of circulation. As will be discussed in detail later, this particu- 
lar state was reached after the vortices were subjected to one or more 
of the demise mechanisms. The fundamental conclusion to be deduced from 
Fig. 23 is that the effect of the stable linear stratification is to 
inhibit the rise of the trailing vortices. This may be due to a number 
of reasons such as: (i) the increase of the downward buoyant force on 

the recirculation cell, (ii) the generation of countersign vorticity out- 
side the recirculation cell (to be discussed later in Section IV. F), 

(iii) the enhanced effect of the entrainment and/or detrainment of the 
recirculation cell, (iv) the possible enhancement of the vortex breakdown 
and alterations in the circulation distribution after vortex breakdown 
due to stratification, (v) the change of the turbulence structure and 
hence the eddy diffusion by stratification, and (vi) the enhancement of 
the Crow instability, or other as yet unknown demise mechanisms associated 
with stratification. 

It is evident from Fig. 23 that the vortices in a stratified medium 
rise to smaller heights. Thus, the free surface effects become important 
at D/b^ values smaller than that corresponding to the unstratified case. 
Clearly, the particular value of D/b^ at which the free-surface effects 
become important in stratified medium depends on the stratification 
parameter. 

The results obtained with Delta 2 are shown in Figs. 23 through 26 for 
three SP values (SP = 0.375, 0.50, and 0.75) and for D/b^ = 6.50. 
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Figure 24„ Vertical Vortex Rise for Delta 2, SP = 0^375, D/b = 6o50„ 
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Figure 25. Vertical Vortex Rise for Delta 2, SP = 0.50, D/b = 6.50 
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Figure 26. Vertical Vortex Rise for Delta 2 , SP = 0.75, D/b = 6.50 



Figure 25 shows the repeatability of the data and all three figures show 
the general consistency of the experiments. 

Figure 27 is a composite plot of the aforementioned figures. A 
perusal of this figure reveals the same conclusions arrived at earlier 
in connection with the discussion of Delta 1 (Fig. 23). 

Since the purpose of the use of two Delta-wing models was to confirm 
the independence of the data from the size of the model, Figs. 28 through 
30 were prepared for comparison of the data obtained from Delta 1 and 
Delta 2. In general, the data compare extremely well for a given strati- 
fication parameter. The minor differences in H/b^ for large values of 
Vot/bo (see Figs. 29 and 30) are thought to be due to the somewhat random 
influence of the demise mechanisms at the later stages of the life of 
the trailing vortices. Finally, it should be noted that D/b^ does not 
affect the variation of H/b with V t/b for D/b larger than at least 
4.3, as discussed earlier. 

D. DEMISE OF DELTA-WING GENERATED TRAILING VORTICES 

The discussion of the demise of the trailing vortices will be carried 
out first with regard to the initial stages of motion and then to the 
fully developed trailing vortices. 

1 . Horseshoe Vortices and Vortex Rings 

As the lifting surface is set in motion impulsively from rest the 
vorticity shed from the surface forms a so-called Horseshoe vortex. 
Shortly after the formation of this vortex, linking occurs and the 
Horseshoe vortex transforms into a vortex ring (see Fig. 31, where the 
Horseshoe vortex is at the point of linking). The observations and 
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Figure 27„ Stratification Effects on Vortex Rise for Delta 



DELTH 1 vs. DELTA 2 SP = 0.375 
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Figure 28. Comparison of the Vortex Rise for Delta 1 and Delta 2 for SP = 0.375. 
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Figure 29o Comparison of the Vortex Rise for Delta 1 and Delta 2 for SP = 0,50^ 
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Figure 30. Comparison of the Vortex Rise for Delta 1 and Delta 2 for SP = 0.75 
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Figure 31. Horseshoe Vortex and the Formation of a Vortex Ring. 





measurements have shown that the vortex ring inclines itself at an angle 
of about 45° and undergoes complex deformations partly under the influ- 
ence of self-induced velocities and partly due to the velocities created 
by the remainder of the trailing vortex. The complexity of the fluid 
motion in the vortex ring is such that analyses based on circular vortex 
rings of constant core and pitch diameter cannot be used to predict the 
transport velocity of the vortex ring. In fact, attempts made to do so 
have resulted in an arbitrary selection of the core radius with widely 
differing velocities. The non-uniformity of the velocities induced on the 
ring resulted in vortex stretching in addition to gross deformations. 

Following the formation of the first ring, the end of the trailing 
vortices necessarily links and forms a new Horseshoe vortex because of 
the fact that the vortices cannot end in a fluid. This, in 
turn, results in a new vortex ring with similar consequences. Repeated 
experiments have shown that as many as three or four such rings are formed. 
Several facts must be brought out in connection with their formation, 
motion, and subsequent demise. Each ring has a character of its own be- 
cause of a number of reasons. The first ring has the largest circulation 
because an accelerating body gives rise to the largest accumulation of 
vorticity. Second, the first ring is formed from that portion of the 
trailing vortex which is relatively young or not too much affected by 
aging. Lastly, the first ring is under the influence of its own self- 
induced velocity plus that of the Horseshoe vortex at the end of the 
trailing vortex. 

The second ring is formed from that part of the vortex trail which 
is relatively more aged. In other words, the aging and deformation of 
a vortex ring formed from a relatively young vortex trail (higher 
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concentration of vorticity at the core) are not the same as that of a 
ring formed from an already aged core. Furthermore, the second ring is 
influenced not only by the velocities induced by itself, but also by the 
velocities induced by the first ring and the remainder of the trailing 
vortex. Extending this discussion to three or four rings one realizes 
that each ring is different at birth and is subjected to different velocity 
fields subsequently. Also, one may surmise that the lack of continuation 
of the formation of a larger number of vortex rings may be related to the 
demise of the trailing end of the vortices by simple eddy diffusion of 
vorticity. 

From the foregoing, one may postulate that the rings will form 
at such a frequency that the total length of the trailing vortex between 
the body and the last ring will remain either constant, or will get larger 
or smaller. Experiments have shown that not only is the number of rings 
finite, but also the rate of horizontal progression. The distance as 
well as the time interval between the successive linkings have been deter- 
mine and it was found that (Ax) /(At) is nearly constant (Ax is the 

n n n 

horizontal distance between the n-th and (n-l)-th link and At^ is the 
corresponding time difference). Comparing the linking velocity given by 
(Ax)^/(At)j^ with the velocity of the model, it has been found that 

(ix)J(At)^ ■ 

Thus, one may conclude that the length of the vortex trail increases 
rapidly and its demise is governed by other mechanisms unrelated to the 
formation of the Horseshoe vortices and vortex rings. However, the 
formation of the vortex rings is extremely important since the formation 
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of the rings is closely related to the change of circulation generated by 
the lifting body. Such a change may result either from the change of the 
angle of attack of the lifting surface and/or from the change of velocity 
of the body. One may also include in the change of angle of attack the 
alteration of the positions of the flaps of an aircraft. Consequently, 
an aircraft or a submerged body moving with rapidly changing velocities 
or at a constant velocity with periodically oscillating lifting surfaces 
can give rise to vortex rings. Finally, it is worth noting that the 
changes in velocity and/or the angle of attack of lifting surfaces may 
precipitate and enhance other demise mechanisms. 

2. Demise of the Vortex Pair 

The maximum height to which the vortex pair rises has been obtained 
from the plotted as well as the tabulated data and is shown in Fig. 32 as 
a function of the stratification parameter Nb^/V^. Only those points for 
which D/b^ were sufficiently large so as to preclude any free surface 
effects are shown in Fig. 32. Also shown in this figure are the couple 
of data points obtained by Tomassian [Ref. 47]. 

Figure 32 points out several important facts. First, the vortices 
have a sel f- limiting effect on their rise, even in an unstratified medium. 
Second, the stratification imposes further limitation on the maximum 
height of rise. In all cases, various demise mechanisms bring about the 
said limitation. These will now be discussed in some detail, bearing in 
mind the fact that it is in some cases extremely difficult to quantify 
some of these mechanisms. 

The most obvious of all the mechanisms is the so-called aging of 
the vortices or the diffusion of vorticity by viscous and turbulent actions. 
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MRXIMUM VERTICRL RISE OF VORTICES 
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Figure 32. Experimental Data for H /b versus Nb /V 



It IS a well-known fact that a vortex with a given initial vorticity 
distribution diffuses its vorticity gradually. If the vortices are 
moving in a stratified medium the density gradient gives rise to opposite- 
signed vorticity outside the recirculation cell. This vorticity may be 
entrained into the recirculation cell, thereby reducing the total circu- 
lation. The entrainment and detrainment of the recirculation cell have 
been discussed in detail by Crow [Ref. 51] and by Tombach [Ref. 59], 
among others, and will not be discussed here further. Suffice it to 
note that the end result of the aging of vortices is to decrease their 
upward migration. 

The next demise mechanism is that first rigorously analyzed by 
Crow [Ref. 54] and known either as Crow instability or sinusoidal insta- 
bility. Crow has shown that two initially parallel trailing vortices are 
unstable to small perturbations and links at intervals of 8.6 times the 
initial core spacing b^. The rings so formed gradually dissipate at a 
rate dependent on atmospheric conditions. In fact. Crow and Bate [Ref. 

55] have shown the role played by turbulence in the sinusoidal instability, 
the wavelength, and the demise of the entire trailing vortex. In spite of 
the uncertainties associated with the quantification of the effects of 
turbulence, the instability analysis of Crow still remains as the only 
rigorous mathematical model of any one of the demise mechanisms. 

In the data shown in Fig. 32, the Crow instability was observed 
in the trailing vortices for stratification parameters less than about 
0.75. At higher stratification parameters, the vortices do not last long 
enough to experience Crow instability. This raises the important ques- 
tion as to the time required from the inception to the linking of the 
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vortices in stratified fluids. Such an analysis does not yet exist and 
the said time will probably be a strong function of the ambient turbu- 
lence as in the case of unstratified fluids. 

The events following the linking of the vortices are both inter- 
esting and complex. The stretching of the vortices, the instability 
propagating about the rings, the increase of the core size, and the 
generation of opposite-signed vorticity give rise to as yet incalculable 
consequences. This aspect of the investigation, like many other inter- 
esting problems raised during the course of this investigation, needs 
further examination, probably through the use of high-speed motion 
pictures . 

The last and certainly one of the most important demise mechanisms 
to be discussed is the vortex breakdown. When the axial velocity in a 
vortex becomes nearly equal to the maximum tangential velocity at the edge 
of the core, the vortex undergoes a rapid transition^ from supercritical 
to subcritical flow. This abrupt transition may occur in one of three 
forms depending on the circulation and the Reynolds number of the ambient 
flow. These are axi symmetric vortex breakdown, spiral breakdown, and 
double-helix breakdown [Ref. 7]. The occurrence of the vortex breakdown 
may be precipitated by an adverse pressure gradient. For example, the 
swirling flow in a conical diffuser produces vortex breakdown at a dis- 
tance must closer to upstream than a uniform pipe. Thus in the case of 
an axi symmetric streamlined body, the deceleration of the ambient flow at 
the aft end produces vortex breakdown in the trailing vortex. In the case 
of the Delta-wings, the breakdown results from the expansion of the core 
and the reduction of the tangential velocity. The consequence of the 
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breakdown is not the total annihilation of the vortex, as it was once 
thought. Observations and measurements have shown that the vortices sur- 
vive the vortex breakdown and continue to rise. It is not yet known as 
to how much circulation if left in the vortex following the vortex break- 
down or a cascade of breakdowns. In all cases, the trailing vortex may 
contain one or more velocity maximums with numerous smaller swirls. It 
appears that the reduction of circulation is more with vortex breakdown 
than with the other mechanisms. 

The data points shown in Fig. 32, particularly for small values 
of the stratification parameter have resulted from the vortices which had 
sufficient time to undergo Crow instability and vortex breakdown. Depend- 
ing on the particular experiment, either the Crow instability or the vortex 
breakdown occurred first since there was no strong adverse pressure gradi- 
ent to precipitate the vortex breakdown first. Since the vortex motion 
resulting from a combination of Crow instability followed by vortex break- 
down is not the same as that resulting from vortex breakdown followed by 
some form of Crow instability, one does expect considerable scatter in 
Hmax/I^o’ evidenced by the data points in Fig. 32, particularly 

for = 0. It must be emphasized that is a measure of the 

maximum relative height at which the vortex cores have essentially dissi- 
pated. It is quite possible that some residual vorticity may still per- 
sist in the flow. This vorticity is not expected to have significant 
consequences relative to the vorticity associated with the trailing 
vortices prior to the dissipation of the cores. In other words, the 
vorticity density in the flow field for H/b^ larger than is insigni- 

ficant relative to H/b^ smaller than 
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The need to determine the lifespan of vortices and to understand 
the physics of the demise mechanisms gave rise to several approximate 
analyses. These will now be discussed in some detail and another 
approximate analysis will be presented. 

Scorer and Davenport [Ref. 50] assumed that the impulse imparted 
to the vortex pair and the circulation of each vortex remain constant. 
This resulted in a gradually decreasing distance between the vortices 
given by; 



V 

H = sinh (ANt) 



( 8 ) 



and 



*^o cosh (ANt) 



(9) 



in which 



A = 0.67 

Eliminating time between Eqs. (8) and (9), one has 

= f no) 

Scorer and Davenport assumed that the recirculation cell will cease to 
rise when the vortices annihilate each other by having their cores touch. 
In other words, assuming b = 2r^, one has from Eq. (10), 
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( 11 ) 



Vx _ 1.492 ^r^oT2 ^ 

b ■ Nb /V \'2r ’ ' ‘ 

The above equation yields the maximum relative height in terms of the 
stratification parameter provided that one has some idea about the core 
radius, or bQ/i^Q- The use of r^ = O.lb^ in Eq. (11) results in unrealis- 
tically large H /b. values. Various values of b /r^ were examined so 
as to obtain some correspondence, at least for some ranges of Nb^/V^. 

This resulted in a choice of b^/r^ = 3. Then Eq. (11) reduces to 

|_| 

max _ 1.67 

^ ■ ^V'^o ^ ^ 

Equation (12) is plotted in Fig. 33 together with the experimental data. 
Clearly, there are two fundamental difficulties associated with Scorer 
and Davenport's analysis. First, the assumed core radius of b^/3 is too 
large to be realistic. Second, as the stratification parameter approaches 
zero, the analysis predicts increasingly large H /b. values, implying 
that a trailing vortex in an unstratified infinite fluid medium will never 
cease to rise. These shortcomings are simple consequences of the fact 
that the physics associated with the demise mechanisms have not been in- 
corporated into the analysis. Judging from the results of other works 
since the pioneering work of Scorer and Davenport, one must state in all 
fairness that it has not yet been possible to incorporate the consequences 
of turbulence. Crow instability, and vortex breakdown into the analysis. 
Attempts made to model numerically to include the effect of turbulence 
will be discussed later. 
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MAXIMUM VERTICAL RISE OF VORTICES 
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Figure 33o Comparison of the Experimental Data with Eqs„ (12) and (15)„ 



Another approximate analysis of the ultimate vertical migration 
of the recirculation cell may be carried out by assuming that the 
initial kinetic energy of the fluid within the cell is converted to 
potential energy. 

The kinetic energy of the fluid within the recirculation cell 
is given by [Ref. 62] 



E 



1 ^ p2,, 0 ^ Iv 

^pr (In — + j) 



(13) 



The potential energy of the recirculation cell may be written as 
2 

in which mb^ is the area of the recirculation cell and m = 2.83. Equat- 
ing the kinetic and potential energies using r^/b^ = 0.1, and simplifying, 
one has 



' iNFTh > 

0 0 

Equation (15) is shown in Fig. 33 together with Eq. (12) and experi- 
mental data. Clearly, Eq. (15) represents the upper bound of the 
H^ax/'^o since no energy dissipation due to any one of the demise 

mechanisms has been incorporated into the analysis. Finally, it should 
be emphasized that the predictions of Eq. (15) could have been partially 
matched to the experimental data by artificially increasing the relative 
core radius r /b^. In the foregoing, was assumed to have the 
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experimentally observed value of 0.1 in order to accentuate the energy 
loss brought about by the complex demise mechanisms. 

A third approximate analysis may be carried out by assuming that 
the vortex separation remains constant, but the strength of the vortices 
vary with time. 

Starting with the analysis of Scorer and Davenport [Ref. 50], 
i.e., by equating the time rate of change of impulse imparted to the 
vortex pair to the buoyant force, one has 

- gmb^^p (16) 



Assuming 3b^/3t = 0, one has 

pb„ ^ ■ g>"b2pBz (17) 

in which iip = pBz and B = N^/g. Noting that ^ ^ = 2 t^’ 

0 

may be reduced to 



rdr 

ub 

0 



g m b„ B z dz 
^ 0 



(18) 



Integration of Eq. (18) and rearrangement in terms of the stratification 
parameter yields 



U 

max 




1.49 

Nb /V 
0 0 



(19) 



The correspondence between the experimental data and the prediction of 
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Eq. (19) (see Fig. 34) may be somewhat fortuitous because of a number of 
reasons. The size of the vortex core does not enter into the analysis 
even though it is known to have an important effect on the life of 
vortices. Furthermore, as in all other analyses, the known demise mechan- 
isms do not appear in Eq. (19). 

All of the foregoing approximate analyses point out emphatically 
that there are no simple approximations with which one can predict the 
rise of the vortices. One must use experimental data and heuristic 
reasoning to delineate the role played by various demise mechanisms in 
the ascent of vortices. This fact is further accentuated by Fig. 35 
where the experimental data are shown together with a best fit curve 
defined by 



= 4.64 exp[-1.3(Nb^/V^)^] (20) 

0 

The particular dependence of on Nb^/V^ strongly suggests that the 

phenomena such as turbulence, entrainment and detrainment. Crow instability, 
and the vortex breakdown play far more complex roles than what may be 
quantified in a simple approximate analysis. 

E, SAILPLANE GENERATED VORTICES 

The results obtained with a streamlined body and its lifting surface 
in stratified and unstratified water will be presented in the following. 

Figure 36 shows the normalized vertical distance H/b^ as a function 
of the relative distance Z’/b^ along the direction of motion for three 
representative model velocities. The reason for the choice of the particu- 
lar normalizing parameters is that the body proximity does not allow 
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Figure 34o Comparison of the Experimental Data with Eq^ (19) „ 
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Figure 35. Experimental Data and the Best-Fit Curve [Eq. (20)]. 
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Figure 36„ The Rise of the Sailplane Generated Vortices as a Function of Z/b 
for SP = 0„ 



the determination of an initial rise velocity V . In spite of the 
scatter, the data show that there is very little or no dip in the path 
of the trailing vortices as they approach the aft end of the body at 
Z/bo = 14. Once the vortices are free from the body proximity effects, 
they continue to rise at a nearly constant velocity. In fact it is this 
particular rise velocity that is used in Fig. 37 to plot H/b^ as a function 
of V^t/b^. 

Also shown in Fig. 36 is the analytical result based on ideal line 
vortices (see Appendix E). The vortex characteristics used in the calcu- 
lations were those appropriate for the lifting surface of the streamlined 
body. The predictions of the inviscid flow model come surprisingly close 
to the experimental data. This may be due to a number of reasons, the most 
important of which is the use in the calculations of the vortex strength 
which corresponded to that experimentally determined value after a series 
of vortex core bursts, i.e., the use of the vortex velocity based on the 
data shown in Fig. 36. 

The exploratory results obtained with one stratification parameter 
(SP = 0.625) are shown in Figs. 38 and 39. A comparison of these figures 
with Figs. 36 and 37 show that the stratification inhibits the rise of 
the vortices as anticipated on the basis of the experiments with Delta- 
wings. The effects of other stratifications on the evolution of trailing 
vortices shed from the lifting surfaces of an axisymmetric slender body 
will have to be examined in light of additional experiments. The determina- 
tion of the circulation existing in the trailing vortices after the vortex 
breakdown will enable one to compare with greater confidence the prediction 
of the numerical models with the results obtained experimentally. 
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Figure 37„ The Rise of the Sailplane Generated Vortices as a Function of V^t/b 
for SP = 0. 
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Figure 38. The Rise of the Sailplane Generated Vortices as a Function of Z/b 
for SP = 0.625. 
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Figure 39„ The Rise of the Sailplane Generated Vortices as a Function of 
V^t/b„ for SP = 0„625„ 



F. NUMERICAL MODEL OF LAMINAR VORTEX PAIR 

As discussed in the previous sections, the need to determine the 
lifespan of vortices and to understand the physics of the demise mechan- 
isms gave rise to several approximate analyses. Numerical methods have 
been employed only recently using computers to attempt to incorporate the 
effect of turbulence in the evolution and ultimate demise of the vortex 
pair. It is hoped that such methods will provide not only a power of 
prediction of the migration of trailing vortices but also some insight as 
to the relationship between turbulence and the demise mechanisms. To be 
sure, there are, at present, no numerical methods dealing with the three- 
dimensional trailing vortices and the effects of vortex breakdown and 
Crow instability. 

The most advanced numerical methods have been applied to the study 
of the two-dimensional, unsteady vortex pair. Hecht et al . [Ref. 44], 
introduced such a model which solves the mean and ensemble-averaged 
Reynolds stress equations of fluid motion, based on a second-order closure 
turbulence model. The initial distribution of vorticity was assumed to 
be Gaussian of the form 

4 = ^ exp[-(r/a)^] (21) 

TTa 

with a = 0.25b (this corresponds to r /b^ = 0.177). The initial turbu- 
lence distribution was again Gaussian with a spread of 0.25b^, and was 
assumed to be isotropic with a representative integral scale. An upward 
velocity was employed to follow the vortex descent. A Reynolds number 

of V b /v = 2x10^ was assumed for the analysis. 

0 0 
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The model was validated by comparing the calculated and experimentally 



determined trajectories of vortex rings . The results of this model for 
a vortex pair in atmosphere with a stratification parameter of 0.8 are 
compared in Fig. 40 with those of the Delta-wings in water with a strati- 
fication parameter of 0.75. The model matched the experimental data very 
well, although there are no physical arguments to justify the selected 
integral scale for turbulence. Furthermore, the value of a should be 
about half that assumed. 

In a recent paper, Hecht, Bilanin, and Hirsh [Ref. 53] compared the 
predictions of the aforementioned model with other analytical [Refs. 49, 

51, and 59] and experimental results [Ref. 43]. Initial distribution of 
vorticity was based on the data from Burnham et al. [Ref. 43]. The initial 
a used in Eq. (21) was determined to be 0.074b^. The integral scale for 
turbulence was set to a smaller value for this analysis and a sensitivity 
study of the macroscale effects on vortex descent was presented. The 
results of this model for SP = 0.67 and a/b^ = 0.075 are compared with the 
Delta-wing results for SP = 0.625 and SP = 0.75 in Fig. 41. This model 
does not match the data for V^t/b^ greater than about 2. 

A better match between the measured and calculated values could have 
been achieved with a = 0.25b^. But this value is larger than that in 
trailing vortices by a factor of about 3 relative to the observations of 
Burnham et al . [Ref. 43]. Of course the integral scale of turbulence 
could have been adjusted also to achieve more realistic results for 
V t/b larger than about 1.5. 

While the numerical models cited above have pioneered in the study of 
turbulent effects in trailing vortices, the difficulties associated with 
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Comparison with Hecht et al Model C443 
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Figure 40„ Comparison of Experimental Data with that Predicted by Hecht et al (Ref^ 44) » 
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Figure 41 „ Comparison of Experimental Data with that Predicted by Hecht et al (Ref. 53). 



the choice of free parameters (a and the integral scale of turbulence) 
point out the need for additional work. 

In an effort to eventually produce such a fully turbulent flow model 
to compare the three-dimensional trailing vortex data, gathered from the 
experimental models in the water basin, a two-dimensional, laminar, un- 
steady flow numerical model of the vortex pair was developed as part of 
this investigation. 

For this purpose, the Navier-Stokes equations were written as 



3u 

3t 
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3t 
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+ vA^v 
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p 



( 22 ) 

(23) 



where the y-axis is assumed to be directed downwards. Eliminating the 
pressure between the two equations, one has 

(24) 

r = ly. _ 3v (25) 

and represents the vorticity. For flows in which the gravitational 
acceleration is several orders of magnitude larger than the fluid acceler- 
ations, the terms in the bracket in Eq. (24) may be neglected. This 
amount to invoking the Boussinesq hypothesis. Thus, the equation of 
motion and the equation of continuity may be written as 
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( 26 ) 



+ iLyil + ^ jy sl = vv^r - -2- ^ 
at ax ay ^ ax 

l£ + iLypl + i(j^ - vv^p ( 27 ) 

at ax ay ^ ^ ^ 

The last term in Eq. (26) represents the effect of the density gradient 
and gives rise to opposite-signed vorticity. 

The equations that were chosen to describe the motion include a 
Gaussian distribution of the initial vorticity given by 



? = ^ exp[-(r^/2r2)] (28) 

Zirr 

0 

Boundary values of the u and v velocity components are determined by 
the Biot-Savart Integrals 

u(x,y) = f ^ — y)^(x »y ) (jx'dy' (29) 

v(x,y) = f (x-x');(x',y') 

•^A Zirr'^ 



Where: 

= (x-x‘)2 + (y-y‘)^ 

x',y' denote the position of vorticity, and 

x,y the point where the velocity components are calculated. 

The stream- function values around the boundaries and velocity within the 
field boundaries are calculated from u and v using 
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V 



dip 

dx 



(31) 



Field St ream- function values are calculated using Poisson's equation 



The calculations had to be carried out in a finite numerical mesh, as 
in all other numerical analyses. The appropriate boundary conditions are 
shown on the boundaries of the quadrant in which the calculations have 
been performed (see Fig. 42). It should be noted that the conditions 
u(0,y) = v(x,0) = 0 are satisfied automatically by evaluating the Biot- 
Savart equations in all four quadrants. Furthermore, the symmetry pro- 
vided by the free surface (x-axis) and by the normal bisecting the vortex 
pair (y-axis) enable one to confine the domain of calculation to a 
single quadrant. 

Additional details of the analysis are as follows: 

(i) The domain of calculation is prescribed; 

(ii) The vorticity distribution given by Eq. (28) is assigned to 
mesh points; 

(iii) A particular stratification is assigned to the density distribution 

(iv) Boundary values of velocities are calculated using the Biot- 
Savart equations; 

(v) The boundary values of the stream- function are calculated; 

(vi) The flow field values of the stream- function are calculated 
using Poisson's equation and a successive overrelaxation scheme; 

(vii) The velocity field is calculated through the use of the stream- 




(32) 



function values. 
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Figure 42. Computational Field and Boundary Conditions for 
the Numerical Model . 
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The foregoing marks the end of the first time step. The successive 
steps are sequenced as follows: 

(i) The vorticity and density fields are calculated using an 
upwind-differencing scheme; 

(ii) The new center of vorticity is calculated using the square of 
vorticity as the weighing parameter; 

(iii) A vertical velocity is imposed on the field using the velocity 
calculated from the rate of change of the center of vorticity; 
and finally, 

(iv) The new points at the top of the field are initialized and the 
vorticity is used to calculate the boundary velocities. 

The last step completes the iteration and marks the beginning of a new 
sequence.. 

The plots of the streamlines of the flow field for a simulation of 
the motion of the vortices generated by Delta 1 for SP = 0 are shown in 
Figs. 43 through 47. The time indicated on each plot is the normalized 
time t/(b^/V^) which is of course identical to its common form V^t/b^. 

In these calculations, D/b^ was taken to be 4.31, and the initial value 
of r^ was assumed to be 0,09b^. The calculated values of H/b^ are com- 
pared with those obtained experimentally in Fig. 48. The matching of 
the two results is surprising in view of the fact that neither turbulence 
nor any other demise mechanism in conjunction with turbulence was taken 
into consideration. The surprising agreement between the two results 
may be attributed to numerical diffusion and to the relative stability 
of the trailing vortices in unstratified fluids. 

Sample streamline plots resulting from the numerical calculation for 
SP = 0.50 are shown in Figs. 49 through 53 for various values of H/b^. 
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Figure 43. Streamlines of the Numerical Model for SP = 0, 

H/b^ = 0. 
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Figure 44. Streamlines of the Numerical Model for 
SP = 0, H/bg = 0.60. 
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Figure 45. Streamlines of the Numerical Model for 
SP = 0, H/b^ = 1.13. 
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Figure 46. Streamlines of the Numerical Model for 



SP = 0, H/b^ = 1.59. 
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Figure 47. Streamlines of the Numerical Model for 
SP = 0, H/b^ = 2.01. 



no 



(O 



in 




O 

-Q 

O 

> 



cvi 



C9 



oq/H 



m 



Figure 48. Comparison of Experimental Data with the Predictions of the Numerical 
Model for SP = 0. 
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Figure 49. Streamlines of the Numerical Model for 
SP = 0.50, H/b = 0. 
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Figure 50. Streamlines of the Numerical Model for 
SP = 0.50, H/b^ = 0.61. 
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Figure 51. Streamlines of the Numerical Model for 
SP = 0.50, H/b^ = 1.17. 
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SP = 0.50, H/b^ = 1.70. 



115 



Y/bo 



SP = 0.50 H/bo = 2.22 



3.82“ A 0PR 

-AA BPBR 
- AA 30HE 

AAA e»3BB 

4 ♦ 4 ^ 0 0 6 

4444 +4+-A 

♦ 4 lUllUl H0 

♦♦ mil nil +A as 

♦ 111 111 ♦ 00 

♦ 11 ? 11 +APH 



♦ 1 y.7222222222 1 1 > 90 

>11 ’>22222222227222 1 > 0 

>11 2222 2222 11>A0D 



ccccccccccrccccccccccrr.c 
ccccccccccccccccccccccc 
ccccccccccccccccccccc 
cccrr^ccci-ccrrcccr.cc 
c cc cc c r c rc c r c ccc c c 
ccccccccccccccccc 
ccccccccccccccc 
cccccc cccccccc 
cccccccrrxccc 
ccc cccccccc 



cccccccccc 
ccccccrcc 
cccccccc 



>1 


222 


222 1 - B 


ccccccc 


>1 


22 


333"»333333 322 11 ♦ 9 


cccccc 


4l 


221 


: 333333 


3333333 22 1> 3 


CCCCC 


1 


22 


■j -3 *3 -a -a "1 


3-»3 -a-iT 3 2’ 1> G 


ccccc 


1 


1 2 


333333^3 


33333333333 222 If 8 


CCCC 


1 


^12 


3333333 


33333333 22 1> 0 


CCC 


1 


22 


333333 


3333333 22 1 8 


CC 


1 


1 2 


333 33 


333333 22 1> P 


cc 


1 


2 2 


33 3 3 3 


33333 222 1 9 


cc 


1 


22 


3333 


33333 22 1 > 9 


c 


1 


2 2 


3333 


33333 22 1 AB 


c 


1 


22 


3333 


33333 22 B 


c 


1 


2 2 


3333 


33333 22 1>B 




1 


1 1 


3333 


33333 22 1> B 




1 


1 1 


3333 


33333 22 U B 




1 


22 


33333 


33333 22 1> 




1 


2 2 


33333 


33333 22 1 -6 




1 


22 ^ 


333333 


333333 22 I A3 




1 


2 7 


-a -a -a 


3333333 722 1 > 0 




1 


22 


33333333 


3*’3 33333 22 1 > 0 


c 


1 


22 


332333333 ’333333333 2 22 1>A8 


c 


>l 


222 


3333333 


33?'’? 323333 222 1 > R 


c 


>1 


222 


3333333 3 '*3 33333 222 llv 0 


c 



>1 

>1 

>11 

>11 
> 1 
> 11 
> 11 
> 111 



222 

222 



333333-2333 
3 33 3 33 



222 

2222 



2722 2222 11 

222222 2^2222 11 
221227272172227222 1 3 

21112222172172 II >A69 

2222222 11 > 9B 

11 <*A 3 



1>A8 0 
1 8 



> 111 111 >A 38 

> nil nil 36 

» iminiiiiiii BBB 

» 111 mill +> A RGB 

»> 1 44 4A BR2 

4444 4444 A BRR 

iA e?r: 

►-AA P3B 
AAA RRB 

— AAA eeeB 
- AAA 38B8 
_ , , - AA BR88 

3.1! — esoee 



3 



3.5 



X/bo 



cc 
cc 
ccc 
ccc 
cccc 
cccc 
ccccc 
cccccc 
cccccc 
ccccc cc 
cccccccc 
ccccccccc 
ccccccccccc 
cccc cccccccc 
ccccc cccccccc 
ccccccccccccccc 
cccccccc cccccccc 
cccccccccc cccccccc 
ccccccccccccccccccc 
cccccccc ccc ccccccccc 
ccccccccccccccccccccc 



T = 2.73A 



Figure 53. Streamlines of the Numerical Model for 
SP = 0.50, H/b^ = 2.22. 
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The comparison with Figs, 43 through 47 show that in the stratified case 
the vortices do not split out laterally as much as in the unstratified 
case. The reason for this is that the opposite-signed vorticity is 
generated, as shown by numerical calculations, outside the lower right 
and lower left sides of the recirculation cell. The effect of this 
vorticity is to pull the vortices together and thus increase the velocity 
of their upward migration. Had turbulent diffusion been present, the 
reduction in circulation would have more than offset the effect of opposite- 
signed vorticity, thereby slowing the rise of vortices. In laminar 
vortices, however, the strength of the vortices changes very slowly and 
the vortices accelerate upward in an amount closely related to the magni- 
tude of stratification. The foregoing arguments, however qualitative, 
point out once again the conflicting conclusions arrived at by various 
workers (see Table I). 

The experimental and calculated values are compared in Figs. 54 and 
55 for SP = 0.50 and SP = 1.0, respectively. The correspondence between 
the measured and calculated values is surprisingly good, partly because 
of the fact that even the experimental values of H/b^ for the stratified 
(SP = 0.50) and unstratified cases do not significantly differ (see Figs. 

23 and 27). However, in the case of SP = 1.0 where the experimental data 
for SP = 0 and SP = 1.0 show significant differences (see Figs. 23 and 27), 
the predictions of the numerical calculation differ significantly from 
the experimental data as seen from Fig. 55. As noted earlier, the strong 
opposite-signed vorticity generated in this case accelerates the vortices 
upward and yields larger H/b^ values. It is interesting to note that at 
the early stages of the motion, where the effect of counter-signed 
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Figure 54„ Comparison of Experimental Data with the Predictions of the 
Model for SP = 0.50. 
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Figure 55. Comparison of Experimental Data with the Predictions of the Numerical 
Model for SP = 1.0. 



vorticity is not yet evident, the slope of the numerical prediction is 
nearly identical to that of the experimental data. 

It should be emphasized that the foregoing was only an exploratory 
numerical investigation of the evolution of a laminar vortex pair. It 
served to confirm the need for the development of a comprehensive 
numerical model which is not only true to the underlying physics of the 
phenomena but also devoid of floating parameters whose values are arbitrarily 
chosen (e.g., initial vorticity diffusion, mesh size, size of 

calculation domain, etc.). The results also show that turbulence plays 
an extremely important role in the diffusion of vorticity and thus in 
the entire life of vortices from their inception to their ultimate demise. 

Even when such a comprehensive model is developed it does not now 
appear that it would be too easy to surmount the difficulties associated 
with the modeling of the consequences of the Crow instability, vortex 
breakdown, and the entrainment and detrainment of the recirculation cell. 

Of course it is not too surprising that such an investigation should 
require a better understanding of turbulence itself through measurements 
and calculations of vortical flows. 

G, SCALE EFFECTS 

The results obtained with Delta wings have shown that the Reynolds 
number does not play a measurable role in the evolution of the turbulent 
trailing vortices. Thus, the data presented herein may be used for larger 
Delta wings provided that the background turbulence is negligible and the 
stratification is linear. Additional experiments are needed to delineate 
the effects of the background turbulence and the nonlinearity of the 
stratification. Such experiments are currently underway. 
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As to the streamlined body, the Reynolds number of the model is about 
800 times smaller than that of a full scale vehicle. In spite of the 
rather large difference between the Reynolds numbers of the model and the 
prototype, however, the data presented herein are thought to be valid for 
the full scale vehicle provided that the background turbulence is negli- 
gible and the stratification is nearly linear. This assertion is based 
on the fact that the boundary layers on the model and its lifting surfaces 
were rendered turbulent with tripping wires and sand strips. Obviously, 
this does not preclude the need for full scale experiments for a number 
of reasons, the most important ones being the effect of background turbu- 
lence and the nonlinearity of the stratification. 

One must bare in mind that there are fundamental differences between 
the demise mechanisms of trailing vortices generated by a Delta wing or 
rectangular wing (in the absence of any body-proximity effects) and those 
generated by a vehicle and its control surfaces. The demise of the vortices 
emanating from a Delta wing is primarily due to sinusoidal instability and 
turbulent diffusion. Thus, the lifespan of these vortices may be strongly 
affected by the background turbulence [Ref, 55], The demise of the trailing 
vortices generated by the lifting surfaces of an axi symmetric body is 
primarily due to the vortex breakdown. In fact, the data presented herein 
show that the vortex breakdown (brought about by the body-proximity effects) 
severely limits the rise of the trailing vortices. The comparison of the 
data obtained with Delta wings and the streamlined body show that the maximum 
rise of the cores of the vortices generated by the sail planes of the stream- 
lined body is about 2,5 times smaller than that of the Delta wings. Thus, 
in the case of vehicles one would be interested not only with the effect of 
background turbulence and the type of stratification but also with the 
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effect of Reynolds number on the intensity and generation of vortex 
breakdown or cascade of vortex breakdowns. The previous studies [Refs, 

6-8] have shown that the occurrence of vortex breakdown is not materially 
affected by the Reynolds number (and Mach number) provided that the 
vortex is initially turbulent. In other words, the background turbulence 
is not of major importance in the occurrence of vortex breakdown, at 
least when the breakdown takes place on a single vortex. One cannot 
ignore the fact that events such as vortex breakdown do not occur simul- 
taneously and symmetrically on a pair of trailing vortices. The mutual 
induction of vortices and the non-symmetrical occurrence of cascades of 
vortex breakdown on each vortex may enhance the role played by the back- 
ground turbulence. Furthermore, the amount of circulation remaining in 
the recirculation cell following the breakdown may be affected by the 
background turbulence. The best estimate, at present, is that the inten- 
sity of the adverse pressure gradient (i,e,, the mean flow field about the 
body) is fundamentally responsible for the occurrence and the position of 
the vortex breakdown. The subsequent demise of the vortices may be affected 
by the background turbulence, mutual induction, and the type of stratifi- 
cation. 

It is evident from the foregoing that it would be very desriable to 
conduct a series of field experiments. In doing so, the data presented 
herein should serve well for their planning and for the interpretation of 
the environmental effects on the evolution of the trailing vortices and 
vortex rings. Clearly, meaningful full scale experiments require, for 
their interpretation, not only the measurement of the characteristics of 
vortices but also the quantification of the prevailing environmental 
conditions. 
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In summary, the evolution of vortices shed by a body and its control 
surfaces in a natural environment with non-isotropic background turbulence, 
currents, shear, wave motion, and stratification is such a complex time- 
dependent phenomenon which demands quantitative analysis, ingenious 
experiments, and qualitative descriptions for the purpose of explaining 
the underlying causes of the physical effects and providing a model by 
which design and operation can be safely and reasonably performed. 
Ultimately, the flow field and the vortex motion about full scale vehicles 
will not be quantified solely due to the knowledge gained through the study 
of well-posed problems and controlled laboratory experiments alone but 
rather through the interaction and reasoble marriage of the two approaches: 
basic studies and full scale experiments. 
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V. CONCLUSIONS 



The investigation reported herein warranted the following conclusions: 

lo In unstratified water, the trailing vortices generated by a 
lifting surface moving at a negative angle of attack rise with time and 
dissipate ultimately at a finite height (about H/b^ = 4.6), The top of 
the recirculation cell rises to a normalized height of about H/b^ = 6, 

This height should not be confused with the height to which the dye in 
the trailing vortex diffuses after a long period of time, 

2, The rate of rise of vortices is always smaller than that predicted 
theoretically through the use of the inviscid line vortices, 

3, The rise of the vortices was shown to be governed by the dimension- 
less parameters V^T/b^, D/b^, and Nb^^/V^ (the stratification parameter) for 
linearly stratified fluids with no background turbulence, 

4, The effect of stratification is to reduce the rate of rise of 
vortices and the maximum height attained by them. In fact, the larger the 
stratification (i,e,, Nb^/V^), the smaller the maximum height attained by 
the vortices, 

5, The ultimate demise of the vortices generated by Delta wings is 
observed to be brought about by sinusoidal instability and turbulent 
diffusion. The demise of the vortices generated by the lifting surfaces 
of an axisymnetric body is brought about by a cascade of vortex breakdowns 
and the turbulent diffusion. 

6, The vortices subjected to vortex breakdown rise to a maximum 
height which is about 2,5 times smaller than those subjected to sinusoidal 
instability, 

7, The results have shown that the body-proximity effects on trailing 
vortices are very strong and always give rise to vortex breakdown. 
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8, Various approximate analytical models were considered and all 
were deemed to be inadequate. None of the approximate models yielded a 
finite rise height in unstratified water and none considered the effects 
of the demise mechanisms, 

9, An exploratory numerical analysis based on two-dimensional, laminar 
unsteady flow assumptions was carried out. The predictions of the model 
were in good agreement in the unstratified case. In the stratified case, 
the generation of the counter-signed vorticity decreased the vortex span 
and increased the velocity of the upward migration of the vortices. The 
predictions of the model could be made more realistic through the inclusion 
of at least the effects of turbulence and possibly those of other demise 
mechani sms, 

10, The scale effects were examined in as much detail as possible. 

It is concluded that the scale effects under similar conditions (linear 
stratification and no background turbulence) are insignificant. However, 
the demise of the vortices undergoing sinusoidal instability may be strongly 
affected by the background turbulence. The body-proximi ty effects and the 
adverse pressure gradient, rather than the background turbulence, are thought 
to govern the occurrence and the position of the vortex breakdowns. However, 
the ultimate demise of the vortices subjected to vortex breakdown may be 
affected by the background turbulence, in addition to other laboratory and 
environmental factors such as non-linear stratification, currents, waves, 
wave-current interactions, wakes of other bodies, shear, etc. 
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APPENDIX A 



RELATIONSHIP BETWEEN THE AMOUNT OF SALT 
AND STRATIFICATION PARAMETER 



The stratification parameter is defined as and is given by 




b I 

_o ^ 

V ^D p 
0 0 



(A-1) 



where: 

N = the Brunt- Vaisala frequency and is defined 'by 

N = <''-2) 

b = the vortex core separation 
0 

2 

g = 32.174 ft/s is the gravitational acceleration 

d = the layer thickness in inches 

= the fresh water density 

AP = the change in density between two layers 

V = the initial vortex rise velocity. 

0 

The ratio of the initial vortex-rise velocity to the model velocity was 
found to be constant for each model, independent of model velocity, i.e., 
V^/U = Constant for a given model. 

The brine-water density p^ in the n-th layer is related to the 
concentration C^ by 
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P„ = P^d t 0.71C„) 



(A-3) 



where the salt concentration C is defined in general by the following 
equation as a function of the amount of salt added; 



1 



W 



1 + 



Fr. water 
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sal t 



(A-4) 



Thus , one has 



or 



4P = P^ . P^_, . 0.71 P„(C„-C„.,) 



- '^n-ll 

0 



(A-5) 



and the Eq. (A-1) becomes. 



Nb„ b 

0 _ 0 

V„ ' (V„/U)U 






71 xgx (C^- 



(A-6) 



Solving Eq. (A-6) with respect to (C^-C^_.j) one has 



, (V /U)U 2 Nb 2 
C_-C_ 1 = ^oTTTg) ^ ^ ^ ^ ^ 



'n n-1 



V 



(A-7) 



where n = 1 ,2,3, . . . . 

Because of the fact that the first layer normally is fresh water 
without salt, i.e., = 0, one has 
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= ( 



H (V /U)U p Nb^ . 

o:f^) E-V-^' (ir)' 



Hence 



^ (V /U)U p Nb^ . 

[-^3^ <ir)^ 

0 0 



C^ = 2C, 



= "'^1 



Also from Eq. (A-4) one obtains 



(A-8) 
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n,salt 
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(A-9) 



Combining Eqs. (A-7), (A-8), and (A-9), one has 



W 



W 



Fr. water 



n,salt 



1 

d EVUlin Nb^T 

d )x[^^^ ]^ x (^)2 



h( 



0.71 g 



(A-10) 



0 0 

Equation (A-10) gives the amount of salt that has to be added to the 
n-th layer for a given value of the stratification parameter Nb^/V^ and 
a given layer thickness d x water weight of fresh water that 

is included in layer of thickness D in the water basin, i.e.. 



’^Fr. water = (36)(3)(d)(62.4) = 561.6 dibs. 
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where d is as previously described. Therefore 



W 



n,salt 



561.6 X d 
1 



nd 



t l2xo:7l g lC-V-] ( 

^ r\ 






- 1 



For example, 

0 0 

'^1 ,salt 
Also, for 

0 0 
^1 ,salt 
^2, salt 



for 

= 0.5, d = 1", U = 2.05 ft/sec, and V^/U 

= 0.01519 lbs of salt. 

= 1.0, d = 1", U = 2.05 ft/sec, and 

= 0.0607 lbs of salt, and 
= 0.1215 lbs of salt. 



(A-11) 



= 0.0391 



= 0.0391 
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APPENDIX B 



TEMPERATURE STRATIFICATION EFFECTS 



A nonuniform temperature stratification which might result from the 
temperature difference between the ambient air and the water in the basin 
could have an adverse effect on the salt stratification and on the 
accuracy of the data obtained. To this end, an analysis has been carried 
out by assuming that the basin is full with water at 50 degree temperature, 
the basin fills at a much faster rate than the rate of heat diffusion, 
there is no heat transfer between the bottom and side walls of the basin, 
and the air temperature at the air-water interface is 65 degrees. 

The governing equation is the one-dimensional unsteady conduction 
equation, given by 



The initial and the boundary conditions are 



T(x,0) = 50° F = Tq 



(B-2) 



T(L,t) = 65° F = T 



(B-3) 
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3T(0,t) 
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(B-4) 



Letting 9 = T - T , one has 

CO 




(B-4) 
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and 
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9(L,t) = 0 


(B-6) 
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A straightforward application of the technique of separation of 
variables yields 



T(x,t) 



T + 2 9^ 
00 0 



n=0 



ci.X^t 
t n 






cos X X 
n 



(B-8) 



where 



(2n+l )tt 
2L 



(B-9) 



For the case under consideration, one has 



t = 4 hrs, L = 45 inches, and = 0.00022 in^/sec 



The numerical results have shown that the temperature distribution in 
the basin is affected to a depth of only seven inches after a period of 
4 hours. The small temperature gradient between the said depth and the 
free surface gives rise to a negligible density gradient in all but the 
regions very close to the free surface, a region which is not of special 
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interest. It is of importance to note that the numerical values used 
in the calculations correspond to those encountered in coldest months 
of the year. For the majority of the test runs the water temperature 
was about 63 degrees compared to the air temperature of 67 degrees. Thus, 
the effect of heat conduction on density stratification was several orders 
of magnitude smaller than that calculated. 

As to the effect of side walls, a lateral temperature gradient may 
develop along the side wall since the thermal conductivity of aluminum 
is much larger than that of water. However, such a temperature gradient 
will have very little or no effect on the vertical gradient and hence on 
the motion of trailing vortices. Experiments with dye particles falling 
slowly in a vertical plane, extending from wall to wall, have not re- 
vealed any natural convection along the side walls. Thus, it is concluded 
that the heat transfer in any direction in the basin did not play an 
important enough role to affect the results. 
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APPENDIX C - DELTA 1 TABULATED DATA 
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APPENDIX D - DELTA 2 TABULATED DATA 
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APPENDIX E 



INTERACTION OF TIP AND ROOT VORTICES FOR 
VARIOUS BODY CONFIGURATIONS 



It is a well-known fact that the combination of the body and lifting 
surfaces may give rise to root vortices in addition to the tip vortices. 
The occurrence of the tip vortices depend strongly on the manner the 
lifting surfaces are attached to the body. If and when they occur, they 
affect the evolution of the tip vortices and hence their ultimate demise. 
Consequently, it is of special importance to examine in detail, albeit 
through the use of the inviscid line vortices, the mutual interaction of 
the two types of vortices for various body shapes. 

The analysis described below assumes three different body shapes: 

(i) a streamlined axisymmetric body with lifting surfaces; (ii) a stream- 
lined half body [half of that used in (i)J and the identical lifting 
surfaces; and (iii) only the lifting surfaces protruding above a plane 
surface (representing a body of infinite radius). 

The case (i) was analyzed through the use of image line vortices and 
the appropriate complex velocity potential (see Fig. 56). The analysis 
yields the following complex velocity 



iu+v 




r2/27T r,/27T 





(E-l) 
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Figure 56. Real and Image Vortices for the 
Whole Body Analysis. 
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where u and v are the x and y components of velocity, -r-| and r^, the 
strengths of the tip and root vortices; z-| and the complex displacement 
vector locations of the vortices from the body center; a, the radius of 
the body (dependent on the distance along the body); and A', the change 
in body cross-sectional area per unit longitudinal length of the stream- 
line cylindrical body. 

Writi ng 
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' ) -(u 
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+ AZ 
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(E-4) 



where At is the time step, the vortex positions were calculated as a 
function of time. Clearly, the values superscripted by (n-1) represent 
those previous to those superscripted by n. The value of A' depends on 
the specific shape of the body under consideration. 

The results are shown in Figs. 57 and 58 for the tip and root 
vortices. In these figures, Y/b^ represents the normalized height of 
the vortices as measured from the streamlined body and Z/b^, the normalized 
distance along the body, starting at the tip of the lifting surface. 

Figure 57 shows the vortex positions shortly after their formation while 
they are under the significant influence of the potential flow generated 
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(Near the Lifting Surface). 



WHOLE BODY RNRLYSIS 
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Figure 58„ Positions of Tip and Root Vortices for the Whole Body 
(Away from the Lifting Surface). 



14 



by the body. It should be noted that the tip vortex dips near Z/b^ = 
where the streamlined body terminates. In other words, the potential 
flow about the aft end of the body pulls the tip vortices toward the 
body. Beyond the end of the body, the four vortices are under each 
other's mutual induction only, without body effects. Figure 58 shows 
the long term interaction between the tip and root vortices (of course 
assuming that the vortices retain their inviscid character). Clearly, 
the tip and root vortices rotate about each other. Real vortices are 
not expected to survive such a continuous rotation. 

The case (ii), (the half-body problem), was analyzed in the same 
manner through the use of the vortex system shown in Fig. 59. The 
following complex velocity was used to obtain the velocity of the tip 
vortex in the first quadrant 
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in which every parameter has the previously ascribed meaning to it. 

The results of the half body analysis are presented in Figs. 60 and 
61. The vortex positions are very close to those calculated for the 
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Figure 59. Real and Image Vortices for the 
Hal f Body Analysis . 
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(Near the Lifting Surface) 



HALF BODY ANALYSIS 
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(Away from the Lifting Surface). 



whole body, particularly shortly after the formation of the vortices. 

Thus, one may conclude from a hydrodynamic as well as practical point 
of view that the half body closely models the whole body while providing 
additional experimental flexibilities (e.g., larger depth to body diameter, 
mounting of the test body, etc.). 

The plane-wall analysis, i.e., case(iii), was carried out through the 
use of simple images as shown in Fig. 62. The following complex velocity 
describes the velocity of the tip vortex in the first quadrant: 



iui +v.| 



r^/27T r^/2iT r^/2iT 



Zi+Zi 



z^-z^ ^l“^^l ^r^2 z^+Z2 Z]-Z2 






(E-6) 



Once again, identical values of the parameters were used in order to 
compare the results of the three cases. 

Figures 63 and 64 show the results of the plane-wall analysis. It 
is clear from a comparison of Figs. 63 and 64 with Figs. 57 and 58 that 
the absence of the body gives rise to significant differences in the 
evolution of the vortices. Also, the dipping of the tip vortex at Z/b^ = 14 
is not possible without the streamlined body. Clearly, the half-body 
configuration models much more closely the whole-body behavior than the 
plane-wall model. 

The foregoing analysis is limited by such assumptions as inviscid 
line vortices, instantaneous vortex roll-up, no boundary layer along the 
body, etc. Nevertheless, the hydrodynamic effects of the body are strong 
enough to discourage the use of the plane-wall model in the experiments. 
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Figure 62. Real and Image Vortices for the 
Plane Wall Analysis. 
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Figure 63^ Positions of Tip and Root Vortices For the Plane Wall Analysis 
(Near the Lifting Surface). 
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(Away from the Lifting Surface). 



The prediction of the spiralling of the tip and root vortices from 
the analysis of the whole or half body required the determination of the 
specific conditions under which the two said vortices cannot spiral 
about each other. Such an analysis was carried out only for the simple 
case shown in Fig. 65 where the body effects are ignored. The use of 
the complex function for the two real vortices and their images together 
with the condition of no spiralling yielded 

a = 8^-^- (E-7) 

1 + 36 ^ 

where a = lr^/T2l snd 6 = x^/x^. Clearly, there exists a region delineated 
by certain values of a and B in which the vortices do not rotate about 
each other. The conditions under which the vortices do rotate about each 
other may be exploited for more rapid demise of the vortices. However, 
it remains to be seen as to how and under which conditions the root 
vortices are generated. This information could only come from experiments 
and is very much related to the lift carrying capacity of the lifting 
surfaces and of the body to which they are attached. 
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Figure 65, The Analysis of No-Rotation Condition of the Tip and 
Root Vortices. 
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